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ABSTRACT 

We present 5 to 15 /im Spitzer Infrared Spectrograph (IRS) low resolution 
spectral data of a candidate debris disk around an M4.5 star identified as a likely 
member of the ~40 Myr old cluster NGC 2547. The IRS spectrum shows a silicate 
emission feature, indicating the presence of warm, small, (sub) micron-sized dust 
grains in the disk. Of the fifteen previously known candidate debris disks around 
M-type stars, the one we discuss in this paper is the first to have an observed 
mid-infrared spectrum and is also the first to have measured silicate emission. We 
combined the IRS data with ancillary data (optical, JHK S , and Spitzer InfraRed 
Array Camera and 24 /xm data) to build the spectral energy distribution (SED) 
of the source. Monte Carlo radiation transfer modeling of the SED characterized 
the dust disk as being very flat (/iioo=2 AU) and extending inward within at 
least 0.13 AU of the central star. Our analysis shows that the disk is collisionally 
dominated and is likely a debris disk. 

Subject headings: Galaxy: open clusters and associations: individual: NGC num- 
ber: NGC 2547, infrared: stars, stars: circumstellar matter, stars: individual: 
2MASS 08093547-4913033 
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Introduction 



The Spitzer Space Telescope (IWerner et al.ll2004l ) has provided the astronomical com- 
munity with a wealth of information on sources with circumstellar disks, both previously 
known and recently discovered. The importance of having a more complete statistical sam- 
ple of sources that have disks is tied in closely with our knowledge of the timescales for planet 
formation. For example, knowing when the inner region of the primordial or protoplanetary 
disk begins to clear itself of dust constrains the lower end of this timescale. Several ground- 
based observational nea r-infrared studies g ive evidence pointing to the clearing of the inner 
disk by 3 to 5Myr (e.g. lHaisch et al.ll200ll ). yet t he current paradig m is that planet forma- 
tion occurs on timescales of 10 Myr or longer (e.g. lHillenbrandll2005l ). The primordial disk is 
essentially an accretion disk (formed as the parental molecular cloud collapsed into a star) 
that is still feeding the newly formed star with material. Mainly, magnetos pheric accretion 



is res ponsible for the clearing of most of the gas and dust in the disk (e.g. iMuzerolle et al. 
20011 ) . however, planet formation also helps with the clearing of the gaseous protoplanetary 
disk by the the formation of planetesimals that sweep up much of the disk's material. The 
formation of larger planetary bodies occurs via collisions of smaller planetesimals (lunar- 
mass oligarchs, see iKenyon &: Bromley! 120061 ; I Chambers! 1200 ll . and references therein) and 
this process generates a secondary dusty disk, mostly devoid of gas, known as a debris disk, 
i.e., the debris disk consists mainly of dust that is produced in the collisions between the 
newly formed planetesimals. Since there need to be colliding planetesimals present to re- 
plenish the dust in the disk, detecting a debris disk is a strong signature of an emerging 
planetary system. 

To constrain the duration of the primordial accretion disk phase and the transition from 
primordial to debris disks it is necessary to analyze stars of intermediate age (between 20 
and 40 Myr). The NGC2547 stellar cluster presents itself as an ideal candidate to study 
the evolution of circumstellar disks beca use of its relatively close distance, 361^g 9 pc, and 



age 38.5_ fi \Mvr (jNavlor fc Jeffries! 120061 1 . To investigate the disk population of this cluster, 



Young et al.l (120041 ) acquired Spitzer data (Program ID 58, P.I. George Rieke) o f the cluster 



NGC2547, consisting of imaging in all four InfraRed Array Camera (IRAC) flFazio et al. 
2004 ) bands and all three Multiband Imaging Photo meter for Spi t zer (M IPS) f lRieke et al. 
20041 ) bands. The preliminary results were reported in lYoung et al.l (120041 ). where they found 
that a very small fraction (~ 5%) of cluster members show evidence for level s of infrared 



excess emission likely due to primordial disks. The other significant result from lYoung et al. 



(120041 ) was the discovery of a candidate late-type star (2MASS 08093547-4913033) that has 
a considerable 24 fim excess. It presents no detectable excess emission shortward of 8 /xm, 
while the 24 um excess is reduced compared to CTTS, suggesting it may possess a debris 
disk with an inner hole. This star is of interest because of the rarity of observations of debris 
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Lestrade et al. 


2006; 


Sieeler et al. 


2007; 


Gautier et al. 


2007 



Forbrich et al.ll2008l ). especially for NGC2547's age. Observing a large excess emission, such 



as that detected towards this particul ar star, could be indi c ative of recent planet formation. 



Numerical calculations carried out by iKenyon &; Bromley! (120021 ) showed that formation of 



planetary bodies that are 1000 km in size or larger produces copious amounts of dust, creating 
greater scattering of light from the central star and producing excess emission. 

The aforementione d IRAC and MIPS data were re-analyzed and recently published 
by iGorlova et al.l (120071 ): from their membership list they found that <1% o f sources show 



IRAC excesses, while 30-45% of B-F members present 24 /im excess emission. IGorlova et al. 



( 120071 ) also discovered an additional M dwarf with excess emission indicative of the pres- 
ence of a debr i s disk . To better probe the fraction of debris disks among lower mass stars 
Forbrich et al.l (120081 ) acquired an extremely deep MIPS dataset (Program ID 21127, P.I. 
Charles J. Lada) and discovered nine additional candidate M-dwarf debris disks. Although 
M-type stars are the most numerous of any type, their debris disk statistics are extremely 
scarce due to their faintness. Acquiring more observations of debris disks around M stars is 
therfore very important to improve our understanding of how these disks evolve and planet 
formation occurs. 

In this paper we report our findings regarding the debris disk around the aforem en- 
tioned member of NGC2547, 2MASS 08093547-4913033, found bv lYoung et "all (booi ) to 
possess considerable exce ss emission at 2 4 /im . This source will be referred to henceforth 
as source #23, following iForbrich et al.l (120081 ) 's nomenclature. To determine the source's 
spectral type, we acquired optical spectroscopy. We obtained Spitzer Infrared Spectrograph 
(IRS) observations to study this source in more detail, in particular, to determine whether 
the disk has silicate emission and to help constrain its spectral energy distribution (SED) 
modeling. The shape of silicate features gives insight into dust grain size and possibly also 
its crystallinity if the signal-to-noise ratio is high enough. 

The organization of the paper is as follows: §[2] summarizes the acquisition and reduction 
of the optical and IRS spectroscopic data, and §[3] presents the results obtained. We discuss 
the implications of our results in §HJ and finally, we list our conclusions in §[5J 



2. Observations and data reduction 



2.1. Spitzer InfraRed Spectrograph data 



The spectroscopic data of source # 23 was acquired on the 19th of June of 2006 (IRS 
campaign 32), through GO proposal Program ID 20750 (P.I. Charles J. Lada). The data are 
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comp rised of low resolution (R ~ 60-127) observations taken using the Spitzer IRS fjHouck et al 



20041 ) Short-Low (SL) module, with a total on-source integration time of 3.2 hours. We built 
four Astronomical Observation Requests (AORs), each with 2 cycles of 240 s of ramp dura- 
tion for the second order of SL (SL2: 5.2-8.7/im) and 10 cycles of 240 s of ramp duration for 
the first order of SL (SL1: 7.4-14.5 /im); the AOR keys are 15096064, 15096320, 15096570, 
and 15096832. We used a pointing control reference sensors (PCRS) peak-up target located 
at (a,5)(J2000) = (08 h 09 m 59.22 s , -49° 161 1.0"), with m v = 8.14 magnitudes. Figure Q] shows 
the orientation of the AORs with respect to the background in NGC 2547 during the obser- 
vations. 



We used the CUbe Builder for IRS Spectra Maps (CUBISM) software flSmith et al. 



20071 ) to reduce and extract the spectrum of the source from the four AORs, for each order. 
The reduction began with the subtraction of the off nod position slit from the on-source 
position slit for each order; this step was meant to (i) eliminate instrumental effects such as 
leftover background gradients from the pipeline flat-fielding, (ii) subtract out the zodiacal 
emission, and (iii) remove rogue pixels (created mostly by cosmic rays). Additional bad 
pixels were identified and removed. Three spectra were extracted for each order, namely, on- 
source sp ectrum, and two loc al background spectra on either side of the source. The SMART 



software (IHigdon et al.ll2004l ) was then used to merge the first and second orders, to construct 
an averaged local background spectrum, and also used to build the final on-source spectrum 
with local background subtraction. Figure [2] shows the average local background used (gray) 
and on-source spectrum with and without local background subtraction. Broad polycyclic 
aromatic hydrocarbons (PAHs) bands are present as absorption features in the averaged local 
background spectrum and on-source with no local background subtraction spectrum - these 
are in fact emission bands present in the off-nod observations (Figure [1] shows faint diffuse 
emission in the off-nod slits). As mentioned above, the off- nod observations were subtracted 
from the on-nod observations, turning the PAH emission bands into "absorption" bands. 
The local background subtraction adds back in the PAH flux removed by the on-off slit 
subtraction (see Figure [2]). From 7.5 /im to 14.8 /im the spectrum has a signal-to-noise ratio 
(SNR) of 33, and from 5.2 /im to 7.5 /im the SNR is 22. 



2.2. Optical spectroscopy at Magellan 



Through Spitzer imaging of NGC2547, 1 Young et al.l (120041 ) detected excess emission at 
24 /im from the candidate low-mass star # 23. We have used an optical spectrum of this 
source to measure its spectral type and to constrain the presence of accretion through the 
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strength of the Ha emission line. The spectrum was obtained on 2004 April 26 during multi- 
slit spectroscopic observations of candidate members of NGC 2547 using the Inamori Magel- 
lan Areal Camera and Spectrograph (IMACS) on the Magellan I telescope at Las Campanas 
Observatory. The spectrograph was operated with the 600 1mm" 1 grism, OG570 blocking 
filter, and slit widths of 0.9", which produced a spectral resolution of 2.0 A. The exposure 
time was 60s. After bias subtraction and flat-fielding, the spectrum of the source was ex- 
tracted and calibrated in wavelength with arc lamp data. We corrected for the sensitivity 
functions of the detectors through observations of a spectrophotometric standard star. The 
optical spectrum of the star # 23 is presented in Figure[3j 



Results 



Through a comparison to standard field dwarfs, we have measured a spectral type of 
M4.5 and an equivalent width W\ = 17 A for H a from the optical spectr um. The spectral 
type we determine agrees with that obtained by I Jeffries fc Oliveiral (120051 ). 



The on-source, local background subtracted, IRS spectrum of Figure [2] w as combined 
with existing multiband wavelength photomet ry, from the optical to MIPS 24 fim (jJeffries et al 
2004! ; lYoung et allbooi iForbrich et aDl2008h to build an SED, shown in Figure SJ The final 
spectrum shows silicate emission, however the variable background and the faintness of the 
source makes it very difficult to draw strong conclusions about the shape of this feature. 

We model the SED using Monte Carlo radiation transfer codes develope d for modeling 



and in terpreting spectra and images of dusty accretion disks and envelopes (see lWhitney et al. 



2003bl laL for a description of the codes ). The codes include multiple anisotropic scattering 
and use the iBjorkman fc Wood! (1200 ll ) algorithm for computing radiative equilibrium dust 
temperatures and thermal spectra. 

T he presence of a silicate feat ure in the SED i ndicates that the disk contains small 



grains (IWeingartner fc Draineil200ll ). so we adopt the iKim et al. 



(119941) interstellar me dium 



Whitney et al.ll2003bl fl. for 



(ISM) dust model for the opacity and scattering properties (see 
a graphical re presentation of th e wavelength dependence of this model). The grain size 
distribution of IKim et al.l (119941 ) model follows a power law (slope ~ -3) between 0.02 /im 
and 0.2 /im, and then falls off exponentially. 

For the object under study we adopt an axisymmetric flared disk density structure and 
assume the disk is heated by starlight. The density structure is parameterized by 



p(r, z) = p exp 



1 /*\ 

2 \h 
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where r is the cylindrical radius, z is the vertical scale, a is the radial density exponent, and 
Po is the normalization constant. The disk scale-height is given by h(r) = h wo (r/100 AU)^ 3 , 
where hioo is the disk scale-height at a radius of 100 AU and (3 = a — 1 is the flaring power. 
The disk extends from an inner radius, -R m i n to an outer radius Rd- 

The lack of longer wavelength data makes it difficult to constrain disk parameters, so 
we fix some parameters based on previous models of irradiated disks. We adopt a surface 
density of the form E(r) ~ r" 1 which is ty pical of models of irradiated disks and fix (3 = 1.25 
and a = 2.25 (e.g. D'Alessio et al. fl998l). T he SED is not very sensitive to the disk outer 



radius (e.g., see models in 



Scholz et al 



20061 ). so we set R d = 300 AU. 



For the central star we use a temperature of T* = 3089, a radius of R+ = 0.47 M , a 
l uminosity of L+ = 0.01 7Lm, and the spectrum is for a 3000 K NextGen model atmosphere 



( IHauschildt et al.lll999l ). With these star, disk, and dust parameters fixed, we constructed 
a small grid of SED models varying R m i n , h wo and total dust disk mass, M d . The SED 
is somewhat degenerate in h±oo and M d and we find a good model fit using /iioo = 2AU, 
M d = 2.7(i?rf/300 AU)Mi unar , and disk luminosity of L disk = 0.035L*. The degeneracies 
in these models may be broken with longer wavelength data which is sensitive to the disk 
mass. We see that the silicate feature is fit rather well using the adopted ISM grain model. 
However, this does not rule out the possibility that grain growth has occurred in the disk 
and the silicate feature is produced by a population of small grains in t he upper disk layers 



with larger grains settled towards the midplane (e.g. IScholz et al.l 120071 ). 



The best fit model is shown in Figure |4] and the corresponding parameters are summa- 
rized in Table [IJ What we can say from the modeling is that the disk is considerably lower 
mass and flatter than T Tauri disks. Our models also indicate that if there is an inner disk 
hole it is not very large, -R m i n < 0.13 AU, smaller than that found in other debris disks, e.g. 
AU Mic, typically many AU. However, the nominal model suggests that the hole is larger 
than the dust sublimation radius. 



4. Discussion 



4.1. The nature of the disk 



Although the results derived from SED modeling; described above are consistent with 
a homologously evolved protoplanetary disk (IWood et al.l 12002a ). the disk appears to be 
devoid of gas and non-accreting: the equivalent width of the Ha emission lin e is 17 A, which 
for th e spectral type in question could originate in the stellar chromosphere fjWhite Sz Basri 
2003h . 



-7- 



To better ascertain the nature of this disk, we first consider the relevant timescales for 
debris disk clearing and apply them for our source. There are four possible mechanisms for 
the removal of dust grains f rom within a debris disk (e.g. iDraine fc Salpeterl Il979l ; iKruegel 
20031 ; iGrigorieva et al.ll2007l ): (i) sublimation, (ii) stellar radiation pressure, (iii) stellar wind, 
and (iv) Poynting-Roberston drag. The sublimation temperatures for dust grains depend 
on their composition: amo rphous carbon grains sublimate at 2000 K, whereas silicate grains 
sublimate at 1500 K (e.g. iPreibisch et al.lll993l ). These temperatures are only reached at 
very close proximity to the stai0; sublimation may be an important factor in the production 
of the very innermost regions of a hole, but it is not a relevant grain removal mechanism for 
the debris disk being discussed. The radiation pressure is also an insig nificant grain removal 



mechanism for debris disks around M mass stars (we refer the reader to Plavchan et a. 



for a detailed discussion on this matter). Regarding stellar wind drag. iPlavchan et al 



2005 



fl2005h 



propose it may be a dominant grain removal process, but only for very young and active M 
dwarfs with extremely high mass loss rates (~ 10 3 times higher than solar or 10~ n M yr _1 ), 
which is not the case for source # 23. Using the mass and radius of Table [TJ we calculated 
the escape velocity for this source to be 299 km s -1 . A Gaussian fit to the Ha emission 
line profile gives a full with at half maximum of 114kms _1 , so it seems winds are not very 
energetic for this source. 

The P-R dra g timescale, i. e., the time is takes for a dust grain to spiral inwards into the 
star, is given by (jKruegelll2003l ): 



2^2 



Tp-R 



4:TTp gr ac r 
3l" 



(2) 



where a and p gr are the radius and density of the dust grain, respectively, and c is the speed 
of light in the vacuum. 

Spitzer IRS spectra of the Deep Impact ejecta of comet 9P/Templel revealed a dust 
mixture composed of amorphous and crystalline silicates, amorphous carbon, water ice, and 
other components , with a mean density of the dust grains for this composition of 2.96 gem -3 



( Lisse et al. 2008hl 2 l. A similar dust composition was found in the debris disk of HD 113766 A 
by lLisse et all (120081 ). who observed the source using 5-35 /im Spitzer IRS data. We assume 



1 using the radial dust temperature profile given by equation 28 from lKamp fc Bertoldil (|2000h . we deter- 
mine that the sublimation radius of silicate grains of average size of 0.2 fim around source #23 is 0.0045 AU 
or 2.0 R*. 

2 the densities of amorphous carbon, graphite, olivine, forsterite, and enstatite are 1.81, 2.16, 3.32, 3.27, 
and 3.20 gem -3 , respectively. 
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the dust in the debris disk of source #23 is similar in composition to that of HD 113766 A 
or comet 9P/Templel and we henceforth use the same grain density for our calculations. 
Equation [2] is then re-written solely as a function of a and r: 



T P „ R = 0.12 



a 



AU 



Myr. 



(3) 



A 0.2 /im-sized dust grain at 1 AU would then take 2.4 x 10 4 years to spiral onto our source. 
On the other hand, a 0.2/im-sized grain located at 41 AU would take ~ 40 Myr to spiral 
onto the star, so small grains located beyond this radius could be leftover material from the 
primordial accretion disk if no other dust removal or reprocessing mechanism is at play. 

Yet there is one more important physical process to be considered: grain-grain collisions. 
These collisions either shatter the dust grains into smaller particles (high collision velocities), 
or help clear the disk by triggering grain growth through agglomeration (low collision ve- 
locities); grain-grain collisions therefore re-distribute the grain size population. Grain-grain 



collis ions will also act to reduce the vertical scale height of the disk (see e.g. iBesla fc Wu 
20071 ). The best-fit SED model for our source does indeed predict a "razor-thin" or very 
flat disk, with a scale heigh t of 2AU at 100 AU (see § E1 Tabled]). The grain-grain collision 
timescale, t co u, is given by (IBackman fc Parescdll993l ): 



Tcoll 



7~orb 



(4) 



where r or b is the orbital period and a is the dimensionless fractional surface density ; the 
latter can be approximated by Ldi S k/L+ (jMoro-Martm et al.l 120071 ; iMoro-Martinl 120071 ) and 
Equation H] is written as: 



/ r \3/2 /M \ 1/2 L± 
T ^0,25(-) (^) — yr (5) 

where r is the distance of the grain from the star and M* is the stellar mass. Combining 
Equations [3] and [5] using the disk luminosity and stellar mass cited in Tabled], and a grain size 
of 0.2 /iin, one finds that the collisional timescale is shorter than the P-R drag timescale for 
a disk radius greater than 2.0 x 10~ 7 AU. In other words, the disk is everywhere collisionally 
dominated. This implies that the disk is not a remnant protoplanetary disk but a bona fide 
debris disk since the dust grains have been reprocessed and are "second-generation" . 

For disks with lower relative luminosities, the inner regions may be dominated by P-R 
drag, however, the outer regions are always collision-dominated. For example, a disk with 
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Ldisk = 3.5 x 10~ 5 L+ with average grain sizes of 0.2 /im, would have an inner hole of (at 
least) 0.2 AU radius because in that inner region the P-R drag timescale is shorter than 
the grain-grain collision timescale, i.e., the dust grains would spiral in towards the star at a 
faster rate compared to the rate at which they are produced by grain-grain collisions. 

To date, there are fifteen known debris disks around M stars, namely, AUMic (GJ803, 
HP 197481), M0.5 V e GJ 182, Ml TWA 7, MO. 5 GJ 842.2, and eleven sources in NGC2547 

20081 ). AU Mic is known to have an inner hole of 17 AU and L<n s k/ = 6 x 10~ 4 



( Forbrich et al. 



( iLiu et al.ll2004l ). Equations [3] and [5] could be used to explain an inner hole of ~0.5AU if 
the grain sizes were 0.01 /im, so in AU Mic's particular system another mechanis m may be 



clear i ng the dust (such as s ubstellar or planetary companion, or stellar wind drag, ILiu et al. 



20041 ; iPlavchan et al.ll2005l . respectively). As for GJ 182, i t has 24 iim exce ss and if an inner 
hole exists it has a radi us smaller than 1 AU according to lLiu et al.l (120041 ) . TWA 7 also has 



24 /im excess emission (ILow et al 



2005T) although i t has not yet been poss ible to determine 



if the disk has an inner hole ([Matthews et al.l 120071 ). iLestrade et al.l (120061 ) determined that 
the disk associated with GJ 842.2 has its peak emission at 140 AU although there is tentative 
25 /im emission indicative of warm dust. 

The measured d isk masses of previo usly known debris are reported as 28 Mi unar (GJ 842.2 
Lestrade et ap2006h . 18Mi unnr C TWA 7lMatthews et alJ l2007h. 2.11M /unar fGJ182 lLiu et al. 
2004 ). and 0.89M^ nar (AUMic Liu et al.l 2004 ). The fractional disk l uminosities, Ld i.«k./L* 
vary between 2xl0 ~ 3 CTWA7 ILow et al.l l2005h and 6xl0" 4 fAuMic Em" et alil2004T ). Re- 



garding disk sizes, ILestrade et al.l (120061 ) determined a radius of 300 AU for GJ 842.2 and 
AU Mic's disk has a radius of 2 10 AU. Finally, there is a wide spread in the ages for these 
sourc es, from 200 Myr fGJ842. 2 ILestrade et allbood ) to 12Myr (AU Mic Istrubbe & Chiang 
20061 ). Although the disk around source #23 has an order of magnitude higher fractional 
luminosity than these other debris disks, it appears to be produced by a comparable amount 
of dust. 

The disk scale height we determined from the SED modeling of the debris disk, 2 AU at 
100 AU, is similar to AUMic disk's scale height, 2.5 AU. These scale heights are mu ch smaller 
than those of T Tauri protoplanetary disks. For comparison, IWood et al.l (l2002bl ) modeled 
the SED of the T Tauri star HH30IRS, using the same Monte Carlo radiation transfer code 
as used in the work presented in this paper, and found that the disk scale height was 17AU 
at a disk radius of 100 AU radius. Another T Tauri star whose SED was modeled using the 
same aforementioned code is GMAurigae, that was found to have a scale height of 8AU at 
100 AU radius. 
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4.2. On the silicate feature and the disk mass 



As mentioned above in §[3], the silicate feature in emission indicates that there are small 
grains in the debris disk. Unfortunately, due to the PAH contamination and the faintness 
of the source, we cannot draw any conclusions regarding silicate crystallization. Although 
we cannot draw more quantitative conclusions from our IRS observations, to our knowledge 
these observations correspond to the first silicate emission detection from a debris disk around 
and M-type star. 

To estimate the temperature of silicate grains we assume their temperature is that of 
a blackbody with an SED peaking at 10 /im, i.e., 290 K. The temperature profile obtained 
from the SED modeling (§[3|) indicates that the disk reaches that temperature at a radius of 
~0.15AU (although the silicate emission could also be arising from the upper layers of the 
disk). The snowline, defined as the locatio n where the temperature corresponds to that of 
water ice sublimation (153 K. lHayashilll98ll ). for this disk would be located at ~0.18 AU. The 
habitable zone, anz, defined as the region in the disk where a planet can sustain liquid water 
at its surface, is cent ered roughly on a radius of 0.13 AU for our source (using the relation 
anz = a/ L+/L q from IScalo et al.ll2007l ). To date there are seven known plan et-host M dwarf 
stars, and these have orbits ranging from 0.02 to 2.3 AU (jEndl et al.ll2008l and references 
therein). So, planet formation can indeed occurr within the habitable zone or within the 
snowline for our source. In other words, there could be terrestrial planet formation occuring 
within the habitable zone, however, more data is needed to study this hypothesis in more 
detail. 

Among the M dwarf debris disk candidate s in NGC 2547 , sour ce ^23 is the one that 
presents the biggest excess emission at 24 /im (jForbrich et al.ll2008l ). Such a large excess, 
indicative of a relatively luminous disk, may be the re sult of a stochastic collision of plan- 
etary bodies within the d isk (e.g. iKokubo et all 120061 ) that triggered a collisional cascade 
(IKenyon fc Bromley! 120021 ). Such a cata strophic event hap pened in our own solar system, 
and it led to the formation of the Moon (IKleine et al.ll2002l ). 



The dust disk mass determined from the SED fitting, 1 x 10 _7 M Q (or 2.7Mi unar ), is 
dependent of the distribution of the grain sizes. The fact that we detect silicate emission is 
an indication that small ISM-like dust grains are indeed present in the disk. Unfortunately, 
no longer wavelength data is available and we are unable to further constrain the grain sizes 
used in our model. We can however assume that the grain sizes are smaller than those 
for debris disks around earlier type stars since radiation pressure is not a significant grain 
removal process for M st ars (e.g., the radiative blowout grain size for an AO- type star is 
7.7/im. iBesla &: Wul 120071 ). Furthermore, in § 14.11 we showed that the grain-grain collisional 
timescale for the system is shorter than the P-R drag timescales, and consequently many 
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small grains, possibly submicron-sized, are produced and remain in the disk. We also showed 
that there is no effective grain removal mechanism for small grains in the disk, for the 
exception of grain growth. We believe our mass estimate is therefore a good approximation 
of the true disk mass for the SED fit we obtain. 



5. Summary 

We have acquired Spitzer low resolution spectral data from 5 to 14 /im of a M4.5 
star, 2MASS 08093547-4913033 or source #23, that is a candidate member of NGC2547, a 
~40Myr old cluster. Combined with previously published optical, near-infrared, IRAC, and 
MIPS data, we built the SED of the source and modeled several stellar and disk parameters. 
We briefly summarize the main results of the paper as follows: 

1. We confirm excess emission longwards of 8/xm for the system and report the first 
detection of silicate emission from a debris disk around an M-type star; 

2. The silicate emission indicates that the dust grains seem to be small: micron or sub- 
micron-sized; 

3. SED modeling describes the debris disk as being very flat (scale height of 2AU at 
100 AU radius), and as having a mass of ~ 2.7 lunar assuming a disk radius of 300 AU. 
The SED model suggests that the disk extends to within 0.13 AU of the stellar surface, 
outside the dust sublimation radius. 

4. Comparison of the various timescales involved in removing and producing dust leads 
us to conclude that the disk is dominated by grain-grain collisions throughout its entire 
extent; 

5. A stochastic collision of planetary bodies could be responsible for the higher disk 
luminosity of this source, and formation of terrestrial planets within the habitable 
zone is entirely possible for this system. 

6. We confirm the spectral type of M4.5 for the source and find that its Ha emission is 
consistent with chromospheric activity; the Ha emission line is also narrow enough to 
rule out strong stellar winds. 

We are very grateful to an anonymous referee for a careful review of the paper. We 
are also indebted to James Muzerolle for building the IRS AORs used, and last but not 
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Fig. 1. — Color composite image where red, green, and blue correspond to MIPS 24 /xm, 
IRAC 8/xm and 5.8 /xm data, respectively. The IRS SL AORs used are overplotted in red 
(1st order) and yellow (2nd order). The target source (2MASS 08093547-4913033 or source 
# 23) is indicated by an arrow. 
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Fig. 2. — Comparison of three spectra: the bottom panel shows the averaged local back- 
ground spectrum (gray) and the on-source spectrum with no local background subtraction 
(black), while the top panel shows the on-source with local background subtraction spec- 
trum. The PAH bands (6.25, 7.7, 8.5, 11.25, and 12.8/im) are marked by dotted vertical 
lines. 




Fig. 3. — Optical spectrum of the source # 23. The inset is a zoom-in on the Ha emission 
line. Based on a comparison of these data to spectra of field dwarf standards, we classify 
this st M4.5. The spectrum is normalized at 7500 A. 



1 10 100 

A( ( u,m) 

Fig. 4. — Spectral energy distribution of source #23. The photometric data points are 
marked as filled circles and are comprised of BVI, 2MASS, IRAC, MIPS 24 /im; the IRS SL 
data is represented by the solid gray line. The black curve corresponds to the best fit model, 
whose parameters are given in Table HJ and the dotted curve represents the photosphere. 
The inset shows a zoom-in on the IRS spectrum and the model fit (the spectrum errorbars 
are represented in lighter gray). 
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Table 1. Stellar and disk parameters for source #23. 



Parameter 


Symbol (units) 


Value 


adopted distance 


D (pc) 


361 


stellar temperature 


r*(K) 


3089 


stellar radius 




0.47 


stellar luminosity 


L* (L ) 


0.017 


stellar mass 


M* (M ) 


0.11 


adopted disk outer radius 


Rdisk (AU) 


300 


disk mass 


M disk (M )t 


lxl0~ 7t 


disk inner radius 


Rhole (AU) 


< 0.13 


scale height at 100 AU 


^ioo (AU) 


2 


disk luminosity 


Ldisk (L^) 


0.035 



Hot a disk radius of 300 AU. 
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